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Abstract 

Palladium is an ideal system for understanding the behavior of hydrogen in metals. In Pd, H is located 
both in octahedral sites and in dislocation cores, which act as nanoscale H traps and form Cottrell atmo- 
spheres. Adjacent to a dislocation core, H experiences the largest possible distortion in a-Pd. Ab initio 
density-functional theory computes the potential energy for a hydrogen in an octahedral site in a-Pd and in 
a trap site at the core of a partial of an edge dislocation. The Pd partial dislocation core changes the environ- 
ment for H, distorting the H-Pd bonding which changes the local potential, vibrational spectra, and inelastic 
form factor for an isolated H atom. The decrease in excitation energy is consistent with experiments, and 
the calculations predict distortions to the H wavefunction. 
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Renewable energy requires new methods for the production, storage and transportation of en- 
ergy from the point of production. The potential for hydrogen as an energy storage medium 1 
has renewed interest in the fundamentals of hydrogen in metals — a topic with a long history.- In 
particular, the ease of catalysis of molecular to atomic hydrogen on the surface of palladium has 
motivated the study of atomic hydrogen and hydride formation in Pd. Hydrogen's vibrational ex- 
citations in metals provide an interesting window into H behavior due to its low mass; and a-Pd is 
a useful system to consider as a model system that is simple to prepare. 3 There are various means 
of measuring the vibrational energies, such as conductance spectroscopy, 4 Raman scattering, 5 and 
inelastic neutron scattering. 6 Neutron scattering measured isotopic effects, 7 linewidths and their 
dependence on temperature, 8 the optical band and its high-energy features, 9 and concentration de- 
pendence of the spectra. 10 Heuser et al. recently measured the hydrogen excitation peak at 4K to 
be 59meV for the dilute (0.08at.%) H concentration, while at 300K the peak was 68meV. n The 
low concentration suggests hydrogen trapped at edge-dislocation core interstitial sites. This is 
due to the strong interactions of dislocations and hydrogen, 12 with a binding energy of 0.2eV. 13 
To determine the effect of the dislocation core on hydrogen vibration spectra, we compute the 
anharmonic H potential energy adjacent to a dislocation core from first-principles, and calculate 
excitation energies and wavefunctions. The result is a decrease in excitation energy consistent with 
experiments and distortions in the hydrogen wavefunctions due to strain and symmetry breaking. 

The anharmonic potential for hydrogen in Pd increases vibrational excitation energies and pro- 
duces anisotropy in the inelastic form factor compared with a harmonic potential. 1415 The low 
mass of H produces zero-point motion of 0.1 5A in an octahedral interstitial environment in a-Pd, 
which in turn samples a potential energy surface that is no longer purely harmonic. This increases 
transition energies for H vibration in Pd by ~ 50% relative to the harmonic approximation. 14 The 
vibrational spectrum for H has the degeneracy of an octahedral potential, 16 and excitations which 
are not simple multiples of the first excitation energy. 10 The wavefunctions have cubic symmetry, 
and variations in the inelastic form factor with scattering wavevector show modulations between 
high symmetry directions.-^ 7 - However, when a hydrogen atom occupies an octahedral site near 
a dislocation core, the cubic symmetry is completely broken, and H experiences a distorted envi- 
ronment. 

Predicting the vibrational excitations for a hydrogen atom requires accurate computation of the 
potential energy for hydrogen. Density-functional theory calculations are performed with vas p 18 ' 19 
using a plane-wave basis with the projector augmented- wave (PAW) method, 20 with potentials 
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generated by Kresse. The local-density approximation as parametrized by Perdew and Zunger— 
and a plane-wave kinetic-energy cutoff of 250eV ensures accurate treatment of the potentials. The 
PAW potential for Pd treats the s- and J-states as valence, and the H s-state as valence. The restor- 
ing forces for H in Pd change by only 5% compared with a generalized gradient approximation, or 
including Pd 4p-states in the valence; our choice of the local-density approximation is computa- 
tionally efficient, and gives an a-Pd lattice constant of 3.8528A compared with the experimentally 
measured 3.87 18A. To compute the dynamical matrix for the lattice Green function 23 , and to relax 
H at the octahedral site in a-Pd, we use a 4 x 4 x 4 simple-cubic supercell of 256 atoms, with a 
6x6x6 k-point mesh; while the dislocation geometry with 382 atoms uses a 1 x 1 x 8 k-point 
mesh. The electron states are occupied using a Methfessel-Paxton smearing of 0.25eV. For the H 
octahedral site in a-Pd and the partial dislocation core, atom positions are relaxed using conjugate 
gradient until the forces are less than 5meV/A. 

First-principles calculations using lattice Green function-based flexible boundary conditions 
compute a stress-free edge dislocation in Pd. Flexible boundary conditions embed a dislocation 
in an infinite bulk which responds harmonically to forces.—" 2 ^ The harmonic response is captured 
by the lattice Green function, which gives the displacement field necessary to relax a line-force in 
an infinite harmonic crystal; it is computed directly from the force-constant matrix. 23 The initial 
geometry is a periodic "slab" supercell of 382 fee Pd atoms: an infinite cylinder with free surfaces 
along the dislocation threading line §[112] in a 9a[110] x 11|[111] box. Anisotropic elasticity 
gives the initial displacements for an §[110] edge dislocation centered in the cylinder. Flexible 
boundary conditions use forces from density-functional theory, and relaxes each atom based on 
its location within one of three different regions determined by the distance from the dislocation 
core. 24,25 The 45 atoms around the core start with non-zero forces, and conjugate-gradient mini- 
mization relaxes them while the other atoms are fixed. The neighboring 100 atoms between the 
core and the free surfaces act as a coupling region to a virtual infinite bulk. They start with zero 
forces, but displacing the core atoms induces forces. With the lattice Green function, we can dis- 
place all 382 Pd atoms to relax the forces in the coupling region. This includes the outer 237 
atoms that have non-zero forces due to their proximity to the free surfaces; their forces are ig- 
nored in the calculation. The relaxation cycle continues alternating between conjugate gradient 
near the core and lattice Green function relaxation in the intermediate region, until all forces are 
less that 5meV/A. The final result is the stress-free dislocation core equilibrium geometry: an edge 
dislocation split into two |(211)-type partials separated by a distance of 6.5b (c.f. Fig.Q])- 26 
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Fig. Q] shows the distorted local environment for a hydrogen atom adjacent to a partial dislo- 
cation core. Interstitial atoms favor regions of expansion, and below the slip plane of the partial 
dislocation core has maximum tensile strain and non-volumetric effects, such as differences in the 
bondlengths along the slip-plane normal. Considering the plane (111) halfway between the three 
Pd atoms above and below, the H atom is displaced away from the partial core (above the top three 
atoms). The three Pd atoms above the interstitial site have longer Pd-H bondlengths than the three 
below it. For an idealized edge-dislocation, the slip-plane normal points along the strain gradient 
from compressive to tensile, and just below the slip plane the predominant strain is uniaxial ten- 
sion along the Burgers vector. Changes in the bondlength above and below the interstitial site are 
greatest for the [100] and [010] octahedral axes, which are non-orthogonal to the Burgers vector. 
This effect is such that although there is an overall 5% volume expansion about the interstitial 
site, the [100] and [010] octahedral axes expand by an average of 3%, and the [001] axis contracts 
by 2%. For the octahedral axes non-orthogonal to the Burgers vector, there are expansion and 
contraction of Pd-H bondlengths above and below the interstitial site, while bondlengths along the 
[001] are approximately the same. 

The potential energy for H in an octahedral site in ar-Pd and a partial dislocation core are fit to 
force versus displacement data for the hydrogen atom from its relaxed position. Displacements of 
0.025 A, 0.1 5 A and 0.2A provide data of restoring force on H, where the 0.025 A displacements 
give the harmonic limit and the larger displacements provides information about the anharmonic 
terms. For ar-Pd and ar-Pd with a 5% expansion in volume, all small displacement directions are 
equivalent; we use [100], [110], and [111] directions for 0.2A displacements for a total of four 
displacements. For the partial dislocation core, cubic symmetry is broken: we use all six (100) 
small displacements and six (100), twelve (110), and eight (111) displacements for a total of 
58 displacements. The potential is represented as a fourth order polynomial in displacements x,; 
excluding linear and constant terms, this gives 31 terms to be fit (reducing to 3 in the case of cubic 
symmetry). This representation of the potential reproduces the first-principles forces to less than 
lmeV/A. 27 

Fig. |2] shows that potentials along directions of the shortened bondlengths are steeper than ar- 
Pd, while those along elongated bondlengths are shallower. For an octahedral site in ar-Pd, the 
(100) directions are the stiffest, while the (111) and (110) directions are shallower. In addition, 
the (110) directions connect neighboring octahedral sites and (111) pass through an octahedral 
face. All of the potentials soften with the expansion of the lattice by 5%. The softest of the 
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(lOO)-type directions have components along the Burgers vector, and are above the interstitial 
site, i.e. [100] and [010]. The resultant potential for coordinates along these directions has large 
cubic terms, reflecting bondlength variation above and below the interstitial site. The potential 
for the [001] coordinate is less cubic, consistent with the smaller change in bondlength along this 
direction. The distortions to the potential from the dislocation core geometry is reflected in the 
broken symmetry for excitations and the resulting wavefunctions for the hydrogen atom. 

Table U shows that the predicted transition energies from our potentials; it is an overestimate 
compared to experimentally measured peaks in scattering intensity, but the reduction in transition 
energy matches well. To solve for the energy spectrum of our potential, we use a basis of products 
of Hermite polynomials in the three displacement directions x; with a fundamental length xo. 
The fundamental length is chosen to harmonically match the anharmonic potential: the harmonic 
ground state |000) is the solution to a ID harmonic potential with stiffness (000|V 2 V|000)/3. For 
a-Pd, x Q = 0.215A and for a partial dislocation core jc = 0.210A; note that these are larger 
than the anharmonic ground-state RMS spread. The full basis \mitn 2 m 3 ) contains all m, up to 
mi + m 2 + m 3 < M, for a total of (M + 3)(M + 2)(M + l)/6 basis functions. Energies for the first 
three excited states are converged to 0.3meV with M = 8 (165 basis functions) and 0.03meV with 
M = 16 (969 basis functions). Our overestimation of the transition energies is consistent with 
other density-functional theory calculations for the octahedral site in a-Pd;— it is also independent 
of the exchange-correlation potential and treatment of the ionic cores. The reduction in excitation 
energy of 7meV in the dislocation core corresponds with the experimental measurement of lOmeV. 

Cubic symmetry is broken at the partial dislocation core, which splits the degeneracy of the 
excited state, and gives three principal axes for the ground state wavefunction spread. The eigen- 
vectors of the spatial covariance matrix, X\. = (xiXj) - (x ; )(x 7 >, give the directions of maxi- 
mum and minimum spread. The ground-state eigenvectors are n\ = [0.66, -0.75, -0.06], ri 2 = 
[0.73,0.66,-0.16], andn 3 = [0.16,0.06,0.99]— roughly, [110], [110], and [001]— with RMS val- 
ues of 0.19A, 0.16A, and 0.14A; compared to 0.15A for H in a-Pd. The soft direction n x « [110] 
corresponds with the softest of the nearly symmetric direction in the potential (c.f. Fig. trian- 
gles); the stiffest direction ri 3 « [001] corresponds with the stiffest nearly symmetric direction in 
the potential (c.f. Fig. [2k, diamonds), and n 2 ~ [110] a mutually orthogonal "bulk- like" potential 
(c.f. Fig. [2]d, diamonds). Moreover, the three directions correspond to the shapes of the first three 
excited states in Table HI 

The first three excited states are p-like, and have maximum spreads along directions which cor- 
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respond to ft\, n 2 , and n 3 . The first excited state has a maximum spread along [0.68, -0.73, -0.07], 
the second along [0.71,0.68,-0.17], and the third along [0.15,0.11,0.98]. The magnitudes of 
maximum spread for each of the three states are 0.35A, 0.27A, and 0.23A, compared to a max- 
imum spread of 0.25A for the first-excited state in a-Pd. The expanded volume in the partial 
dislocation core decreases the transition energy; however Table U shows degeneracy-splitting due 
to broken symmetry from the local strain. This strain allows for a low lying transition state ft\ that 
is 30meV below the similar transition in a-Pd, or 20meV below in the expanded lattice. This is a 
direct consequence of the hydrogen occupying a site near the partial dislocation core. 



Fig. |3] shows the orientation of the first three excited states by plotting the inelastic form fac- 
tor, S (q) from the ground state to each excited state. The inelastic form factor, for an excitation 
from the ground state |0) to excited state S(q) = |(0|exp(-z'g • r)|n,)| 2 determines the inten- 
sity of inelastic scattering with scattering direction q. In Fig. [31 we rotate q from n x to n 2 to n 3 
and back to ft\ . For comparison, we also compute S (q) with bulk excited states projected along 
Hi. The excited states have similar structure to the excited states in a-Pd site. The changes in in- 
elastic form factors are due to the expansion of the excited state along the it,- directions compared 
with bulk, where expanded wavefunctions decrease S(q). Density-functional theory calculations 
combined with accurate treatment of boundary conditions compute the relaxed geometry for a 
hydrogen atom adjacent to a partial dislocation core in Pd. We extract the local potential energy 
for the hydrogen atom from first-principles to predict the quantum mechanical transition energies, 
and compare to the transition energies in an octahedral site in a-Pd. The changes in the excitation 
energies can be directly traced to the distortions in geometry around the hydrogen atom, the po- 
tential energy, and ultimately the wavefunctions for hydrogen. The predicted decrease due to the 
volumetric expansion is similar to the change in experimental measurements at low temperatures 
where the hydrogen occupancy of partial dislocation cores should come to dominate. This pro- 
vides an important step in understanding the changes in vibrational spectra from hydrogen in Pd 
as temperature changes and the hydrogen atoms migrate to more energetically favorable sites in 
the sample. In addition, the potential provides a starting point for considering pipe diffusion of H 
along dislocation cores in Pd. 
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TABLE I. Calculated and measured transition energies for H in a-Pd and adjacent to a partial dislocation 
core. Hydrogen occupies an octahedral site in a-Pd; the transition from the ground state to the triply- 
degenerate first excited state is larger than the experimentally measured peak. A volumetric expansion 
of 5% reduces the predicted transition energy by 9meV. The partial dislocation core produces a similar 
volumetric expansion, but includes additional distortion that lift the degeneracy of the first excited state, 
giving a set of three low-lying transition energies. The labels ft, correspond to the different orientations of 
the excited state wavefunction. 





Theory [meV] 


Experiment [meV] 


a-Pd 


87 


69 [10] 
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78 




«i 


56 




ft 2 


80 


59 [VjJ 


*3 


102 
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b=[110] 

FIG. 1. (a) Strain at octahedral sites in Pd edge-dislocation core; (b) Hydrogen (red) in the octahedral site 
below the Pd (blue) partial dislocation core in the (111) slip-plane. The edge dislocation has a total Burgers 
vector |[110], while the partials have Burgers vector |[211] and f [12l] and are separated by 6.5b. (b) The 
hydrogen is in an octahedral site with large tensile strain; the top three Pd atoms are the most expanded, 
and the bottom three are less expanded. The dislocation core displaces the neighboring Pd atoms for the 
octahedral site from the relaxed bulk positions shown by the (100) arrows. Moreover, the hydrogen atom 
moves away from the top three atoms ([100], [001], and [OlO]) increasing the distance to 1.96A-2.16A. 
The hydrogen moves closer to the other three atoms, decreasing the distance to 1.87A-1.91A, compared to 
a distance of 1.96A for H in an octahedral site in fcc-Pd. 
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x[A] x[A] x[A] x[A] 



FIG. 2. Hydrogen potential energy as a function of displacement along (a) (100), (b) (111), (c) (110) 
directions in the slip plane, and (d) (110) directions out of the slip plane. Thick black shows response 
in a-Pd, and thin black shows response in 5% expanded a-Pd. Positive (negative) displacements along 
corresponding directions are dashed (solid). Arrows show corresponding displacements, for the potential 
curves of the same color and linestyle. Thick displacement arrows have positive out-of-page components; 
negative for thin arrows. Single points on the plots are the results of individual first-principles calculations, 
and the continuous curves are fits to the Hellmann-Feynman forces, (a) The cube directions [100] (circles), 
[010] (squares), and [001] (diamonds) show the greatest change to the potential due to the dislocation 
core. The softest directions, [100] and [OlO], point towards regions of expansion, while the [100] and 
[010] stiff directions bring the hydrogen closer to Pd atoms, (b) The octahedral directions [111] (circles), 
[111] (squares), [ill] (diamonds), and [111] (triangles) show less change due to the dislocation core. The 
most significant softening is for the octahedral direction corresponding to the slip plane normal, (c) The 
out-of-plane closed-packed directions [101] (squares), [011] (diamonds), and [1 10] (triangles) are similarly 
softened in the direction of expansion and stiffened in the opposite direction compared with bulk, (d) The 
in-plane closed-packed directions [101] (squares), [110] (diamonds), and [Oil] (triangles) are all stiffened 
compared with bulk. 
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rotation angle (degrees) 



FIG. 3. Inelastic form factor S(q) for a H adjacent to the dislocation core and in octahedral bulk sites. The 
momentum direction q rotates between the three orthogonal vectors corresponding to the directions of 
maximal (it\) and minimal ($3) spread of the H ground state wavef unction in a dislocation core. The peak 
along each n, indicates that the excited state is p-\ike along that direction. The calculation of form factors for 
bulk excited states which maximize overlap with the dislocation core states shows a similar peak shape with 
direction. The height of the peaks corresponds to the spreading of the excited states along those directions 
in comparison with the ground state: increased spreading of the excited states along the particular direction 
decreases the peak height. 
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TABLE S 1 : Pd atomic positions relative to a H atom, and H restor- 
ing force versus displacement at an octahedral site near a partial 
dislocation core in Pd. All distances are in A and forces in eV/A. 



Pd[*] Pdty] Pd[z] Pd[r] 
+2.1562 -0.1321 +0.0016 2.1602 
-1.8553 -0.1860 -0.1090 1.8678 
+0.2754 +1.8880 -0.0351 1.9083 
+0.1341 -2.1461 -0.0170 2.1504 
+0.1869 -0.1148 +1.8790 1.8918 
+0.1788 -0.1726 -1.9491 1.9649 





H [6y] 


H[tf z ] 


H [\6\] 


F x 


F y 


F z 


+0.0250 








0.0250 


-0.0209 


-0.0035 


-0.0028 


-0.0250 








0.0250 


+0.0436 


+0.0115 


+0.0064 





+0.0250 





0.0250 


-0.0009 


-0.0206 


+0.0021 





-0.0250 





0.0250 


+0.0141 


+0.0216 


+0.0010 








+0.0250 


0.0250 


+0.0019 


+0.0045 


-0.0389 








-0.0250 


0.0250 


+0.0112 


+0.0034 


+0.0384 


+0.1500 








0.1500 


-0.0855 


-0.0289 


-0.0189 


-0.1500 








0.1500 


+0.4072 


+0.0722 


+0.0430 





-0.1500 





0.1500 


+0.0353 


+0.0618 


+0.0004 





+0.1500 





0.1500 


-0.0680 


-0.2842 


+0.0089 








+0.1500 


0.1500 


-0.0468 


+0.0245 


-0.3651 








-0.1500 


0.1500 


+0.0127 


+0.0170 


+0.2447 





+0.1061 


+0.1061 


0.1500 


-0.0677 


-0.1353 


-0.1881 





-0.1061 


-0.1061 


0.1500 


+0.0391 


+0.0543 


+0.1692 


-0.1061 





-0.1061 


0.1500 


+0.2519 


+0.0524 


+0.1548 


+0.1061 





+0.1061 


0.1500 


-0.0891 


-0.0091 


-0.2527 


-0.1061 


-0.1061 





0.1500 


+0.2781 


+0.0719 


+0.0274 


+0.1061 


+0.1061 





0.1500 


-0.0993 


-0.2083 


-0.0085 



(continued on next page) 
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TABLE SI (continued) 



H[<y 


H [Sy] 


H[tfJ 


H [|<5|] 


F x 


F y 


F z 





+0.1061 


-0.1061 


0.1500 


-0.0302 


-0.1488 


+0.1392 





-0.1061 


+0.1061 


0.1500 


-0.0032 


+0.0533 


-0.2292 


-0.1061 





+0.1061 


0.1500 


+0.1834 


+0.0542 


-0.1575 


+0.1061 





-0.1061 


0.1500 


-0.0637 


-0.0161 


+0.1592 


-0.1061 


+0.1061 





0.1500 


+0.16259 


-0.07984 


+0.02941 


+0.1061 


-0.1061 





0.1500 


-0.05727 


+0.03962 


-0.01581 


+0.0866 


+0.0866 


+0.0866 


0.1500 


-0.1007 


-0.1365 


-0.1711 


-0.0866 


-0.0866 


-0.0866 


0.1500 


+0.2233 


+0.0673 


+0.1352 


-0.0866 


+0.0866 


+0.0866 


0.1500 


+0.0911 


-0.0486 


-0.1086 


+0.0866 


-0.0866 


-0.0866 


0.1500 


-0.0428 


+0.0367 


+0.1322 


+0.0866 


-0.0866 


+0.0866 


0.1500 


-0.0658 


+0.0377 


-0.2014 


-0.0866 


+0.0866 


-0.0866 


0.1500 


+0.1380 


-0.0549 


+0.1176 


+0.0866 


+0.0866 


-0.0866 


0.1500 


-0.0796 


-0.1479 


+0.1133 


-0.0866 


-0.0866 


+0.0866 


0.1500 


+0.1708 


+0.0664 


-0.1310 


-0.2000 








0.2000 


+0.6714 


+0.1104 


+0.0657 


+0.2000 








0.2000 


-0.1022 


-0.0349 


-0.0233 





+0.2000 





0.2000 


-0.1121 


-0.4854 


+0.0136 





-0.2000 





0.2000 


+0.0388 


+0.0781 


+0.0012 








-0.2000 


0.2000 


+0.0032 


+0.0314 


+0.3769 








+0.2000 


0.2000 


-0.0801 


+0.0408 


-0.5928 





-0.1414 


-0.1414 


0.2000 


+0.0409 


+0.0612 


+0.2424 





+0.1414 


+0.1414 


0.2000 


-0.1040 


-0.2078 


-0.2711 


-0.1414 





-0.1414 


0.2000 


+0.3813 


+0.0782 


+0.2023 


+0.1414 





+0.1414 


0.2000 


-0.1028 


-0.0059 


-0.3836 


-0.1414 


-0.1414 





0.2000 


+0.4267 


+0.0791 


+0.0407 


+0.1414 


+0.1414 





0.2000 


-0.1174 


-0.3327 


-0.0091 





+0.1414 


-0.1414 


0.2000 


-0.0544 


-0.2320 


+0.1817 



(continued on next page) 
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TABLE SI (continued) 



H[<y 


H [Sy] 


H[<y 


H [|<5|] 


F x 


F y 


F z 





-0.1414 


+0.1414 


0.2000 


-0.0173 


+0.0566 


-0.3468 


+0.1414 





-0.1414 


0.2000 


-0.0727 


-0.0158 


+0.2308 


-0.1414 





+0.1414 


0.2000 


+0.2705 


+0.0784 


-0.2222 


+0.1414 


-0.1414 





0.2000 


-0.0678 


+0.0472 


-0.0200 


-0.1414 


+0.1414 





0.2000 


+0.2386 


-0.1188 


+0.0417 


-0.1155 


-0.1155 


-0.1155 


0.2000 


+0.3284 


+0.0769 


+0.1775 


+0.1155 


+0.1155 


+0.1155 


0.2000 


-0.1211 


-0.2055 


-0.2430 


+0.1155 


-0.1155 


-0.1155 


0.2000 


-0.0517 


+0.0435 


+0.1885 


-0.1155 


+0.1155 


+0.1155 


0.2000 


+0.1240 


-0.0680 


-0.1449 


-0.1155 


+0.1155 


-0.1155 


0.2000 


+0.1948 


-0.0769 


+0.1444 


+0.1155 


-0.1155 


+0.1155 


0.2000 


-0.0794 


+0.0425 


-0.2982 


-0.1155 


-0.1155 


+0.1155 


0.2000 


+0.2458 


+0.0737 


-0.1821 


+0.1155 


+0.1155 


-0.1155 


0.2000 


-0.0964 


-0.2252 


+0.1497 
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TABLE S2. Pd atomic positions relative to a H atom, and H restoring force versus displacement at an 
octahedral site in a-Pd. All distances are in A and forces in eV/A. 



Pd W 


Pd[y] 


Pd[z] 


Pd[r] 


+ 1.9555 








1.9555 


-1.9555 








1.9555 





+1.9555 





1.9555 





-1.9555 





1.9555 








+ 1.9555 


1.9555 








-1.9555 


1.9555 



H[<y 


H [Sy] 


H[<y 


H [\6\] 


F x 




F z 


+0.0250 








0.0250 


-0.0239 








+0.2003 








0.2003 


-0.3489 








+0.1418 


+0.1418 





0.2005 


-0.1723 


-0.1723 





+0.1157 


+0.1157 


+0.1157 


0.2003 


-0.1217 


-0.1217 


-0.1217 


+0.2897 


+0.2897 


+0.2897 


0.5018 


-0.3745 


-0.3745 


-0.3745 



TABLE S3. Pd atomic positions relative to a H atom, and H restoring force versus displacement at an 
octahedral site in ey = 5% expanded o--Pd. All distances are in A and forces in eV/A. 

Pd[jc] Pd[y] Pd[z] Pd[r] 
+ 1.9849 1.9849 

-1.9849 1.9849 

+1.9849 1.9849 
-1.9849 1.9849 
+1.9849 1.9849 

-1.9849 1.9849 



HfcJ 


H [Sy] 


H[«y 


H [\6\] 


F x 




F z 


+0.0254 








0.0254 


-0.0129 








+0.2036 








0.2036 


-0.2479 








+0.1441 


+0.1441 





0.2038 


-0.1098 


-0.1098 





+0.1176 


+0.1176 


+0.1176 


0.2036 


-0.0725 


-0.0725 


-0.0725 
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TABLE S4. Polynomial expansion for potential energy for hydrogen fit to force versus displacement data. 



V a -Pd(r) = + 0.5004(x 2 + y 2 + z 2 ) + 4.6069(x 4 + / + z 4 ) - 3.7379(x 2 .y 2 + y 2 z 2 + z 2 x 2 ) 
V +5 %a-Pd(r) - + 0.2487(x 2 + y 2 + z 2 ) + 4.3358(x 4 + / + z 4 ) - 2.2241(x 2 .y 2 + y 2 z 2 + z 2 x 2 ) 
V coie (f) = + 0.6397x 2 + 0.4187.V 2 + 0.7677z 2 + 0.3074x.y - 0.0248.yz + 0.1870zx 

- 2.3703x 3 + 1.6500/ + 0.8837z 3 

+ 0.9825xy 2 + 1.0302z 2 x - 0.8675x 2 .y - 0.8461jz 2 - 0.4979zx 2 - 0.1783.y 2 z 
-0.1212xyz 

+ 4.0708x 4 + 3.5534/ + 5.5368z 4 

- 2.4136x 2 y 2 - 3.2983y 2 z 2 - 3.4471z 2 x 2 

+ 1.6832xy 3 + 0.5148z 3 x + 1.3317x 3 .y + 0.0190jz 3 + 0.8594zx 3 - 0.1617;y 3 z 
+ 0.3463x 2 jz - 0.3546xy 2 z - 0.9699xjz 2 
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